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Abstract: 

We present a new determination of the solar nitrogen abundance making use of 3D hydrodynami- 
cal modelling of the solar photosphere, which is more physically motivated than traditional static ID 
models. We selected suitable atomic spectral lines, relying on equivalent width measurements already 
existing in the literature. For atmospheric modelling we used the CO 5 BOLD 3D radiation hydrody- 
namics code. We investigated the influence of both deviations from local thermodynamic equilibrium 
("non-LTE effects") and photospheric inhomogeneities ("granulation effects") on the resulting abun- 
dance. We also compared several atlases of solar flux and centre-disc intensity presently available. As 
a result of our analysis, the photospheric solar nitrogen abundance is A(/V)=7.86 + 0.12. 

Keywords: hydrodynamics - line:formation - radiative transfer - Sun: abundances - Sun: granulation - 
Sun: atmosphere 



1 Introduction 

Images of the Sun reveal a very complex surface struc- 
ture, which is referred to as granulation, and may be 
understood as the signature of the convective motions 
in the photosphere. Traditional static ID model photo- 
spheres ignore all this complex phenomenology. In the 
last ten years hydrodynamical simulations of stellar pho- 
tospheres have considerably improved and are now at the 
stage that they can be compared realistically against ob- 
servations. This class of models (here and after referred 
to as "3D models") is physically better motivated, al- 
though computationally considerably more demanding, 
than traditional static ID models. For the Sun, the com- 
parison of present 3D models with observations shows 
encouraging agreement (e.g. Figure [TJ. The application 
of 3D models for abundance work is a largely unexplored 
territory, but promising work is in progress not only for 



the Sun but also for other solar-type stars. We began con- 
centrating on chemical abundance determinations based 
on the analysis of high resolution spectra and the use 
of 3D models. The solar abundances clearly occupy a 
prominent place in this project, since we are able to ob- 
tain spectra of very high resolution and S/N ratio for the 
Sun. The pioneering works in this field by Allende Pri- 
eto et al. d200lh and Asplu nd et alj d2004h have led to 
a substantial downward revision of the solar metallicity, 
which implies an a wkward tension with the helioseismic 
measurements (see lBasu & Antiall2008L and references 
therein). It is thus not unsurprising that we start our in- 
vestigation by a reassessment of the abundances of the 
main contributors to the solar metallicity, Z, i.e. oxygen, 
nitrogen and carbon. The present contribution reports on 
a redetermination of the solar nitrogen abundance based 
on atomic lines. 
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Figure 1: Left: Quiet solar granulation as observed with the lm Swedish Solar Telescope (courtesy Mats 
Carlsson 2004). Right: High-resolution C0 5 BOLD simulation of solar surface convection. Both images 
show the emergent continuum intensity (using identical scaling) at A 4364 A in a field measuring 15" x 15" 
(11x11 Mm). 



In Asplund et al. I d2005l) . among other elements, the 
solar abundance of nitrogen was also considered, based 
on a 3D model. Unfortunately not much information was 
given as to which lines are used, oscillator strengths and 
other details of the analysis. Both molecular and atomic 
lines were considered, giving an abundance which is 7.73+ 
0.05 and 7.85+0.08, respectively (including NLTE-corrections) 

2 Models and 

line formation codes 

Our analysis is mainly b ased on a 3D model com puted 
with the CO 5 BOL D code dFrevtag et al.l2002ll2003l: Wede- 
meyer et al. l2004h . Some basic information on the setup 
of this numerical simulation can be found in Caffau et 
al. d2008D . who used the same solar model to determine 
the solar oxygen abundance. We just point out the basic 
differences between the new approach using 3D models 
and the old approach using ID models: while ID mod- 
els describe a time-independent, hydrostatic atmosphere, 
a 3D model is the result of solving numerically the time- 
dependent hydrodynamic equations together with the equa- 
tion of radiation transport. For any instant of time ("snap- 
shot" in 3D-jargon) the 3D models give the physical quan- 
tities on a 3D mesh of points in the photosphere. This al- 



lows a more realistic description of the atmosphere (see 
FigureflJ, since both vertical and horizontal fluctuations 
of physical quantities can be taken into account. More- 
over, the (turbulent) velocity field in the stellar atmo- 
sphere is automatically obtained without the need to spec- 
ify free parameters (like micro- and macroturbulence). 
.We would like to point out that, although not requiring 
the free parameters needed in ID computations to adjust 
the efficiency of the convective energy transport and the 
strength of the turbulent velocity field, 3D models are 
characterized by a set of numerical parameters, e.g. the 
numerical scheme used for solving the hydrodynamical 
equations, the spatial resolution of the numerical grid, the 
amount of artificial viscosity, the number of rays (angles) 
considered in computing the radiation field, and many 
others. The hope is that the results become essentially 
independent of the choice of these parameters once the 
numerical resolution exceeds some critical threshold. 



Besides the CO 5 BOLD model, we considered also 
several ID models for comparison. T hese include the 
semi-empirical Holweger-M iiller model ( Hoi we g ej 19671 : 
iHolweger & Miilledll974 here after HM), a ID model 
computed with the LHD code (see Caffau & Ludwid 20071 
for further details), an ATLAS solar model computed 
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by F. CastellQ and several ID models with solar pa- 
rameters com puted by ourselves with version 9 of the 
ATL AS code jKuruczlll993al l2005bl) in its Linux ver- 
sion dSbordone et al J 120041 : ISbordonell2005l) . using the 
"NEW" Opacity Distribution Functions (ODFs; Castelli 
& Kurucz 2003), and, finally, a ID model obtained by 
temporal and spatial averaging of the 3D model over sur- 
faces of equal (Rosseland) optical depth, which we call 
(3D) model. 

The 3D spectrum synthesis computations are all per- 
formed with Linfor3CQ, which can also compute line for- 
mation using different kinds of ID models as input. For 
comparison, in the case of ID models we also used the 
WIDTH code for calculating equivalent widths (Kurucz 
11993d l2005bl: ICastellil 12005. ISbordond 120051) and the 
SYNTHE code in its Linux version .Kurucz. 1 993b, 2005b; 



ISbordone et al.l2004l : [Sbordonel2005l) for calculating syn- 
thetic spectra. 



3 Line selection, atomic and 
observational data 

In the literature one can find different choices for the 
log gf values (see Table [TJ. We decided in favour of 
the NIST data, first of all because, in this database both 
computed and measured values are critically examined; 
second because all elements ara available, and third to be 
consistent with the other papers we produced on photo- 
spheric solar abundances. 

Our chosen line list and oscillator strengths are given 
in Table [2] 



Results from the line profile fitting 
using different solar atlases 

We also derived the nitrogen abundance from line pro- 
file fitting for few selected lines.We did this exercice for 
comparing the abundance derived from the four high res- 
olution, high S/N, solar atlases available, the two centre- 
disc intensity atlases (the "Delbouille" atlas, i.e. Del- 
bouille et al. 1 19731 and [Delbouille et al.lll98ll and the 
"Neckel" intensity atlas. lNeckel & LabUl 19841) as well as 
the two solar fl ux atlases (the "Kurucz" solar flux atlas, 
Kurucz 2005a and the "Neckel" solar flux atlas, Neckel 
& Labs 1 19841) . From previous investigati ons of the solar 
observed spectra (see iCaffau et al.ll2008l) we know that 
these atlases do not always agree. In Figure [2] we can 



see that the two centre-disk solar atlases are not in agree- 
ment for the line at A 821.6 nm. The solar nitrogen abun- 
dance could be deduced from line profile fitting, but the 
twelve selected nitrogen lines are blended with molecu- 
lar and atomic transitions. We are not sure of the oscil- 
lator strength of the blending lines and we do not have 
NLTE computations for these blending components. But 
for comparing the results from the four solar atlases this 
is not a problem. The four lines we chose are the clean- 
est, although also for these lines blends or close-by lines 
are evident (see Figure[3. . For this reason, being this one 
only a comparative analysis, we performed this exercise 
only with ID models, for which the complete line list 
with blends is available. 

With the ATLAS, HM, and <3D> solar models as in- 
put to SYNTHE, changing the nitrogen abundance, we 
computed three different grids of synthetic spectra. The 
fitting code, described in lCaffau et alj20 05. is based on a 
X 1 minimisation, and for this purpose uses the MINUIT 
procedure. 



1 .00 F 




http://wwwuser.oats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 
http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 
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Figure 2: Observed profile of the Ni line at 
A 82 1 .6 nm as extracted from the Delbouille (solid) 
and Neckel (dashed) solar disc-centre intensity 
spectral atlases. The differences between the two 
solar atlases are not easily explained; possibly tel- 
luric absorption affect the Neckel profile. 



Table [3] shows that the scatter in the abundance de- 
rived from the 82 1 .6 nm line from the four solar atlases is 
considerably larger than what is obtained from the other 
two lines. The systematic uncertainty due to the choice of 
a specific solar atlas that affects the abundance measure- 
ment is computed by comparing the average abundance 
obtained considering each atlas. This uncertainty is on 
average ± 0.02 dex. We did not use 3D synthetic spec- 



4 



Table 2: Selection of Ni lines in the optical and 
near IR bands 



A (nm) 


Ei„„ (eV) 


log gf 


Q 


744.229 


10.330 


-0.385 


B+ 


746.831 


10.336 


-0.190 


B + 


821.634 


10.336 


+0.132 


B + 


822.314 


10.330 


-0.271 


B + 


868.340 


10.330 


+0.087 


B + 


871.883 


10.336 


-0.336 


B + 


1010.513 


11.750 


+0.219 


B + 


1011.248 


11.758 


+0.607 


B + 


1011.464 


11.764 


+0.768 


B + 


1050.700 


11.840 


+0.094 


B 


1052.058 


11.840 


+0.010 


B 


1053.957 


11.844 


+0.503 


B 



Notes:Error of log gf < 0.03 dex (Q=B+), < 0.08 dex 
(Q=B). 



tra for the fitting procedure. In fact, some tens of lines 
would be necessary to consider the whole range. 3D run 
are very time consuming and are, for the time being, not 
able to handle too many lines. Nevertheless, we com- 
pared the 3D synthetic profile to the observed solar spec- 
tra. Two examples are visible in Fig.[3] Only nitrogen 
is considered in the 3D profile; for this reason the syn- 
thetic profile is not able to reproduce the complete shape 
of the feature (see line at 746.8 nm) or gives a too high 
abundance (see line at 868.3 nm). 

4 NLTE computations 

The emergent flux depends on the LTE assumption. This 
is a good approximation for lines forming deeply in the 
photosphere, where the density is high. In fact LTE is 
valid if the photon mean free path is shorter than the dis- 
tance over which matter temperature varies significantly. 
The photon mean free path depends on the probability 
for the photon to be thermalized, hence on the rate of the 
collisions between the absorbing atoms and electrons or 
hydrogen. This rate increases with matter density. 

Since we do not yet have a code able to solve the 
NLTE problem for nitrogen in the case of a 3D model, we 
computed the departures from LTE for the (3D) and the 
HM model using the Kiel code ( Steenbock & Holweger 
1 19841) . with the model atom of iRentzsch-Holml dl996h . 
To take into account excitation and ionisation of the ni- 
trogen atoms by inelastic collisions with neutral hydro- 
gen atoms, the Ki el code uses a g eneralisation of the for- 
malism found in Drawin (1969). This formalism intro- 
duces a scaling factor, Sh, that permits to modify the ef- 
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Table 3: A(N) of Ni from line fitting using three 
different ID models 



Observed 


A 




A(N) L te 




A(N) L te 


Spectrum 


(nm) 


HM 


ATLAS9 


<3D> 


Average 


KF 


746.8 


7.861 


7.824 


7.809 


7.873 


NF 




7.874 


7.837 


7.827 


7.899 


NI 




7.890 


7.834 


7.823 


7.922 


DI 




7.908 


7.849 


7.840 


7.883 


Average 




7.883 


7.836 


7.825 




Scatter 




0.020 


0.010 


0.013 


HM 


KF 


821.6 


7.847 


7.811 


7.805 


7.841 


NF 




7.918 


7.870 


7.867 


7.861 


NI 




7.944 


7.879 


7.873 


7.864 


DI 




7.808 


7.752 


7.746 


7.829 


Average 




7.879 


7.828 


7.823 




Scatter 




0.063 


0.059 


0.060 


ATLAS9 


KF 


868.3 


7.910 


7.889 


7.857 


7.824 


NF 




7.906 


7.877 


7.857 


7.850 


NI 




7.933 


7.878 


7.875 


7.857 


DI 




7.932 


7.886 


7.882 


7.823 


Average 




7.920 


7.883 


7.868 




Scatter 




0.014 


0.006 


0.013 


<3D> 



Notes: Col. (1) is the observed spectrum identification, 
KF: Kurucz Flux, NF: Neckel Flux, NI: Neckel Intensity, 
and DI: Delbouille Intensity; col. (2) is the wavelength; 
col.s (3)-(5) the A(N)lte from line profile fitting with 
HM, ATLAS9, and (3D) model, respectively; col. (6) av- 
erage values of the three lines for each atlas, using the 
model indicated in the last line of each block. 

ficiency of collisions with hydrogen atoms (0 < Sh < 1). 
Currently we do not know which value of the scaling fac- 
tor is the correct one, and therefore decided to compute 
the NLTE-corrections for the two extreme cases (Sh = 
and Sh = 1), and an intermediate value largely used in 
the Kiel community, Sh = 1/3. It could well be that 
each nitrogen transition we considered, actually require 
a different value of Sh- Since, in any case, the NLTE- 
corrections are small, we shall consider this differential 
effect as negligible. 

We considered the effects of the horizontal temper- 
ature fluctuations on the NLTE-correction. The proce- 
dure we used is similar to the one that Aufdenberg et 
al. ( 120050 used to estimate the effects of horizontal 
temperature inhomogeneities. We ordered the emerging 
flux as a function of temperature. We divided in 12 bins 
in increasing temperature and produced horizontally and 
time-averaged models. We computed the NLTE correc- 
tions for these twelve average models. 

From the results of this 1D-NLTE computation for 
each line we found that the NLTE-corrections are small 
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Figure 3: Observed profiles of two Ni lines are 
compared with a 3D synthetic profile. 



and not exceeding -0.05 dex. Regarding the effects of the 
horizontal temperature inhomogeneities, we found that 
lines with lower excitation energy are more sensitive than 
lines with higher excitation energy that are formed deep 
in the photosphere. But these effects are small for our 
sample of nitrogen lines, since all of them have high ex- 
citation potential (cf. TableO. From Figuref4]we can see 
that the horizontal variation of the NLTE corrections is 
0.05 dex at most. For this reason we expect that a full 3D- 
NLTE computation would not differ from our 1D-NLTE 
calculation by more than 0.03 dex. 

5 Nitrogen abundance 

Our final 3D-NLTE nitrogen abundance is obtained by 
averaging the individual 3D-NLTE abundances of each 
line with equal weight. These abundances are obtained 
with the spectrum synthesis code Linfor3D (see Tabled. 
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Figure 4: For two representative lines the NLTE 
corrections from the twelve group-averaged mod- 
els, ordered according to increasing continuum in- 
tensity from left to right, are shown together with 
the result for the global (3D) model (rightmost, 
bold symbol). The lowest curve is plotted at the 
true ordinate level, while the other is shifted up by 
0.20 dex for clarity. 

In the same table we show the total 3D-corrections (de- 
fined as A(N) 3D -A(N) 1Dlhd ) and the granulation correc- 
tion (defined as A(N) 3D -A(N) <3D >) for each line. The lat- 
ter corrections (col. (12)) are negative for all lines, indi- 
cating that the horizontal temperature fluctuations sys- 
tematically strengthen the (high-excitation) lines. The 
total 3D-corrections (col. (13)) are systematically more 
positive, because the IDlhd model produces slightly stronger 
lines than the (3D) model. The 3D-corrections are of the 
same order of magnitude as the NLTE-corrections. 

Our final result for the solar photospheric nitrogen 
abundance is: 

A(N) = 7.85+0.12 for S H = 

A(A0= 7.86 + 0.12 for S H =1/3 (1) 

A(N) = 7.87+0.12 for S B = 1 

If only the EWs from lBiemont et al.l dl990t> are consid- 
ered, A(N) is the same while the scatter is reduced to 
0.06 dex. 

Discussion 

In theory, the same nitrogen abundance should be derived 
from each line. In practice, this is not the case due to 
uncertainties in the analysis, related to the following el- 
ements: the model atmosphere, the log gf values and the 
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EWs. The results discussed in the following refer to the 
3Dnlte model, assuming S H = 1/3. 

Concerning the log g/values and referring to Table 
[2] one could select only lines with the Q value equal to 
B+ and consider the abundances derived only from these 
lines. From the results of Table [4] we find that selecting 
only lines with Q = B+ does not reduce the scatter (both 
for results from the 3D and the HM model) and the mean 
value becomes A(A0=7.83 dex, very close to our recom- 
mended value, A(N)=7. 86 dex. 

' BiemqnTe^alJl|T990) are slightly differ- 
iGrevesse et al.l (1990) by 1 to 6 percent 



The EWs of] 
ent from that of ( 
at most. This is an indication that these EWs are reli- 
able. A possible way to try to decrease the scatter is to 
select the lines for w hich we have the EWs bo th from 
iGrevesse et"al] jl990t) and lBiemont et alj .19901) . Since 
both the two authors chose these lines, we can consider 
this subset more reliable. We find that the mean value 
does not ch ange, while the scatter for the nitrogen abun- 
dance from IGrevesse et al.U 1990) decreases to 0.07 dex. 
A different possibility is to take the four lines, namely 
the lines 821.6nm, 871.8nm, 1011.2nm and 1011.4nm 
for which EWs of the previous authors differ by less than 
three percent. This agreement may be taken to imply that 
the EWs of these lines are more reliable than the others. 
The res ults from this subset are A(A0=7.88 ± 0.06 with 
EWs of lGrevesseeTai](ll990h andA(A r )=7.89±0.06 with 
EWs of Bie mont et 

aLfll990). 

From this subset, it can 
be reasonable to discard the line at 101 1.4 nm. In fact 
we could not obtain a good fit for this line and the cor- 
respondent abundance was too high if compared with the 
fitting results of the other three lines (see Section 3). This 
behaviour is the same for each model and each atlas we 
used. The nitrogen abundance we obtain from these three 
lines is A(N) = 7.86 ± 0.05 dex with EWs of Grevesse 
et al. dl990k and A(N ) = 7.87 ± 0.05 dex with EWs of 
iBiemont et all (}990). With this selection the scatter is 
further decreased for both sets of EWs. 

We conclude from this exercise that the mean nitro- 
gen abundances given above are robust against exclusion 
of the lines which have either the more uncertain log gf- 
values or exclusion of the lines for which we consider the 
EWs to be less reliable. 

Our preferred value for the solar nitrogen abundance 
is A(/V)=7.86±0.12, assuming the NLTE-correction with 
Sh = 1/3. The indicated error represents the line-to- 
line scatter and does not include any uncertainty in the 
log g/values. 
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Table 1: Compilation of the log gf values of permitted N i lines used by the different authors. 



A 


NIST 


Lambert 68 


Lambert 78 


Biemont 


Grevesse 


Takeda 


Rentzsch -Holm 


(nm) 




CA 


NBS 




L 


V 








744.2 


-0.385B+ 


-ft 33 


-ft 4S 


-0.33 


-ft ^87 

U.JO / 


-ft 4(S^ 


-0.411 


-0.386 


-0.573 


746.8 


-0.190B+ 


-0.16 


-0.27 


-0.16 


-ft 171 


-ft 948 


-0.208 


-0.171 


-0.397 


818.4 


-0.286B+ 


—ft 9^ 


—ft AD 














820.0 


-1.001B+ 












-0.996 


-1.017 


-1.090 


821.6 


+0.132B+ 


+ft 1 3 


-ft D1 

\f.\t 1 


+0.13 


+ft 146 


+ft DRQ 


-0.106 


+0.147 


+0.012 


822.3 


-0.271B+ 












-0.288 


-0.267 


-0.390 


824.2 


-0.256B+ 












-0.260 


-0.262 


-0.360 


859.4 


-0.334B 


-0.32 


-0.38 


-0.32 


-0.320 


-0.332 








862.9 


+0.075B 


+0.08 


+0.03 


+0.08 


+0.090 


+0.078 


+0.082 


+0.090 


+0.069 


865.5 


-0.627B 


-0.62 


-0.65 


-0.62 


-0.603 


-0.616 


-0.608 


-0.603 


-0.630 


868.3 


+0.087B+ 


+0.11 


-0.05 


+0.11 


+0.115 


+0.102 


+0.109 


+0.116 


-0.051 


870.3 


-0.321B+ 


-0.29 


-0.41 


-0.29 












871.1 


-0.234B+ 


-0.18 


-0.34 














871.8 


-0.336B+ 


-0.26 


-0.43 


-0.26 


-0.338 


-0.347 


-0.344 


-0.337 


-0.419 


904.5 


+0.439B 












+0.430 


+0.429 




939.2 


+0.320B 


+0.31 


+0.24 


+0.31 


+0.328 


+0.378 


+0.354 


+0.328 


+0.316 


1010.5 


+0.219B+ 












+0.220 


+0.234 


+0.200 


1010.8 


+0.431B+ 


+0.39 


+0.41 


+0.39 


+0.443 


+0.420 


+0.431 


+0.443 


+0.403 


1011.2 


+0.607B+ 


+0.58 


+0.60 


+0.58 


+0.622 


+0.600 


+0.611 


+0.623 


+0.588 


1011.4 


+0.768B+ 


+0.74 


+0.76 


+0.74 


+0.778 


+0.755 


+0.766 


+0.778 


+0.751 


1050.7 


+0.094B 












+0.249 


+0.249 


+0.250 


1052.0 


+0.010B 












-0.045 


+0.010 


-0.040 


1053.9 


+0.503B 


+0.52 


+0.51 








+0.529 


+0.525 


+0.530 


1075.7 


-0.608C+ 












-0.098 


-0.098 


-0.080 


1238.1 


+0.247C+ 












+0.284 


+0.175 


+0.320 


1246.1 


+0.480B 












+0.463 


+0.437 


+0.451 


1246.9 


+0.629B 












+0.622 


+0.622 


+0.610 



Notes: In Lambert 68: CA theoretical Coulomb approximation; NBS (National Bureau of Standards) is the old NIST; in 
Biemont: L length formalism, V velocity formalism. 
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Table 4: A(N)lte, A(N)nlte, and 3D-corrections for selected N i lines with EWs from the literature using 
Linfor3D 



A 


EW 


A(N) LTE 3D 


A(N) NLTE 3D 


A(N) NLTE 3D 


A(N) NLTE 3D 


3D-<3D> 


3D-lD LHD (aMLT=l-0) 


(nm) 


(pm) 






Sh = 


1/3 


Sh 


= 1 


Sh 


= 


dex 


dex 




G B 


G 


B 


G 


B 


G 


B 


G 


B 






(1) 


(2) (3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


744.2 


0.26 0.27 


7.808 


7.826 


7.774 


7.792 


7.782 


7.800 


7.770 


7.788 


-0.039 


-0.002 


746.8 


0.52 0.49 


7.961 


7.931 


7.923 


7.893 


7.931 


7.901 


7.919 


7.889 


-0.033 


+0.008 


821.6 


0.86 0.87 


7.854 


7.860 


7.802 


7.808 


7.817 


7.823 


7.790 


7.796 


-0.039 


+0.007 


822.3 


0.24 


7.593 




7.554 




7.565 




7.544 




-0.055 


-0.018 


868.3 


0.78 0.81 


7.828 


7.849 


7.781 


7.802 


7.794 


7.815 


7.767 


7.788 


-0.037 


+0.007 


871.8 


0.42 0.43 


7.927 


7.939 


7.887 


7.899 


7.898 


7.910 


7.875 


7.887 


-0.044 


-0.006 


1010.5 


0.18 


7.956 




7.939 




7.944 




7.931 




-0.066 


-0.020 


1011.2 


0.35 0.36 


7.897 


7.912 


7.878 


7.893 


7.883 


7.898 


7.869 


7.884 


-0.060 


-0.011 


1011.4 


0.55 0.54 


7.976 


7.966 


7.955 


7.945 


7.961 


7.951 


7.937 


7.927 


-0.053 


-0.001 


1050.7 


0.14 


8.002 




7.992 




7.995 




7.986 




-0.064 


-0.020 


1052.0 


0.08 


7.829 




7.819 




7.822 




7.812 




-0.067 


-0.024 


1053.9 


0.32 


7.989 




7.978 




7.980 




7.971 




-0.057 


-0.010 


average 




7.885 


7.890 


7.857 


7.862 


7.864 


7.871 


7.848 


7.851 


-0.051 


-0.008 


scatter 




0.114 


0.053 


0.122 


0.060 


0.120 


0.058 


0.122 


0.059 







Notes: Col.s with G are from lGrevesse et all jl990h . with B are from lBiemont et all (T990). Col. (1) is the wavelength; 
col.s (2)-(3) the EWs; col.s (4)-(5) A(N) LTE from 3D model; col.s (6)-(ll) A(N) NLTE from 3D model for various values 
of Sh; col.s (12)-(13) two different 3D-corrections. 



